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基于 /2.1 范 数 原子 选择 的 图 像 分 块 稀 琉 重 构 
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摘 要 : 针对 图 像 压 缩 采 样 中 原子 的 选择 规则 难以 确定 的 问题 ， 在 改进 的 正 交 匹配 追踪 算法 的 基础 上 提出 了 一 种 基于 
[1 范 数 的 原子 选择 方式 。12.1 范 数 的 原子 选择 方式 考 卡 了 原子 间 的 相关 性 ， 别 除了 干扰 原子 ,选择 出 了 代表 性 原子 。 将 
所 提 方 法 用 于 图 像 分 块 重 构 ， 算 法 以 图 像 进行 分 块 ， 利 用 P.1 范 数 选择 对 图 像 块 支撑 集 进行 第 选 ， 增 强 块 特征 的 判别 
性 ， 提 高 原子 的 稀疏 度 ， 最 终 提高 图 像 重 构 的 准确 率 和 速率 。 实 验 结果 表明 ， 相 同 条 件 下 在 保证 重建 速度 的 同时 ， 所 
提 新 方法 提高 了 图 像 重 构 精度 。 
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Image segmentation sparse reconstruction based on /2.1-norm stomic selection 


Zhu Hua!’, Yue Jun!', Li Zhenbo’, Zhang Zhiwang!, Kou Guangjiel 
(1. College of Information & Electrical Engineering, Ludong University, Yantai Shandong 264025, China; 2. College of 


Information & Electrical Engineering, China Agricultural University, Beijing 100083, China) 


Abstract: In order to solve the problem of selecting atoms for image compressed samples automatically, a new /2.1-norm atomic 
selection method is proposed. The proposed method takes into account the correlation between the atoms, eliminating the 
interference of atoms and selects a representative of the atomic. Image reconstruction based on image segmentation and [2,1- 
norm selection. In this algorithm, image is divided into blocks, and the image blocks supporting set is filtered by the /2,1-norm 
selection which enhances the discriminant of block features, improves the sparsity of atoms, and ultimately improves the 
accuracy and speed of image reconstruction. The simulation results show that the reconstruction accuracy can be improved with 


the same condition. 
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日 一 化 完备 字典 ， 在 一 定 的 相似 性 度量 基础 上 ， 和 迭代 地 从 字典 
中 选择 最 佳 原子 ， 得 到 稀疏 解 。 在 贪 禁 迭代 算法 中 ， 有 正 交 匹 
前 ， 压 缩 感 知 (CS)〉 叫 作为 信号 重 构 的 新 方法 得 到 极 大 配 追 踪 算 法 [I (orthogonal matching pursuit,OMP) 、 正 则 化 正 
关注 。 因 为 该 理论 打破 了 传统 的 奈 奎 斯 特 抽样 定理 ， 提 高 了 数 交 匹 配 追 踪 贸 (regularized orthogonal matching pursuit , ROMP) 
据 的 传输 效率 。 压 缩 感知 基本 假设 是 : 当 系 统 表示 是 欠 定 时 ， 以 及 采用 回溯 思想 的 压缩 采样 匹配 追踪 算法 外 (CoSaMP) 和 子 
即 观 测 向 量 的 维 数 m 远 小 于 原始 信号 的 维 数 n ， 稀 下 信 号 空间 追踪 算法 0 (SP) ， 以 上 算法 都 需 预测 信号 的 稀疏 度 。 针 
xzeR" 可 以 从 稀疏 测量 矩阵 中 中 重 构 出 。 其 中 ， 原 始 信号 可 能 对 稀 玻 度 未 知 的 情况 ，2008 年 Do 等 人 提出 自 适 应 匹配 追踪 算 
不 稀疏 但 是 在 某 个 变换 域 下 是 稀 琉 的 。 法 0 (sparsity adaptive matching pursuit, SAMP ) ,随后 又 出 现 了 
言 号 重 构思 想 是 通过 最 优化 方法 从 低 维 数据 中 最 大 程度 的 稀疏 度 自 适 应 子 空间 追踪 算法 (sparsity adaptive subspace pursuit, 

和 

林 


并 


恢复 出 高 维 数据 ,压缩 感知 理论 可 使 数据 采集 和 压缩 同时 进行 ， SASP)W23 、 正 则 化 自 适 应 匹配 追踪 〈fregularized adaptive 
对 于 信和 号 的 重 构 有 重要 意义 。 其 中 代表 性 的 重建 算法 有 贪 焚 迭 matching pursuit，RAMP ) 算法 0]、 分 段 正 则 化 正 交 匹 配 追 踪 算 
代 算 法 中、 约束 优化 加、 近似 优 化 外 、 同 伦 算法 外 。 贪 栖 迭 代 算 法 0405、 基 于 约束 等 距 性 的 稀 玻 估计 匹配 追踪 算法 09、 基 

法 所 是 其 中 应 用 较 多 的 一 类 方法 ， 其 核心 思想 是 在 迭代 中 痪 历 标 下 降 和 双 残 差 极 小 化 的 自 适 应 匹配 追踪 算法 上 9。 在 研究 和 总 
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法 的 基础 上 , 本 文 将 4 范 数 和 


压缩 采样 的 思 


本 文 提出 的 方法 可 
1 ”压缩 感知 


假设 有 一 


全 问 量 vA 


以 提高 图 像 的 重 构 精 度 。 


E 论 分 析 与 仿 


图 像 重 构 算法 中 , 提出 了 基于 改进 正 交 匹配 追踪 算法 的 4 
原子 选择 的 图 像 信号 重 构 算 法 。i 


个 实 值 信号 X〔〈 和 eq ) ， 可 以 用 


进行 线性 表示 ,于 


X= 


如 果 X 是 稀 玻 的 , 则 表明 X 是 可 压缩 的 。 如 果 用 


想 融 入 到 
21 范 数 
结果 证 明 


Nxl 维 的 规 


是 任意 信号 X 都 可 以 表 


(1) 


个 与 变 


各 自 表示 。7; 为 算 阵 m 的 i 行 j 列 的 元 素 ，m; 表示 和 矩阵 的 第 


朱 ” 华 ， 等 : 


基于 12,1 范 


行 的 行 向 量 。 
过 2 范 数 0 选择 原子 建立 的 模型 为 

min|> 一 4X|， 本 了 |X|， (8) 

其 中 : +|XD, 为 正则 化 项 ，> 为 平衡 参数 ， 在 分 类 或 拟 合 中 ， 

用 来 避免 过 拟 合 或 欠 拟 合 。 在 此 ， 对 正则 化 项 的 拟 合 系数 进行 
0. 范 数 约束 四。 具体 计算 公式 如 下 : 

[Xb =2 5 -Zl (9) 

该 范式 可 以 在 样本 拟 合 原子 的 时 候 , 通过 对 X 行 向 量 之 间 


换 基 灵 不 相关 的 测量 矩阵 @ (M*N,M<<N) 对 @ 进行 线性 变换 ， ”的 1 范式 约束 M < N 尽 可 能 将 无 用 原子 对 应 的 拟 合 系数 约束 
得 到 测量 向 量 Y: 为 零 。 而 对 列 范式 的 约束 则 保证 所 有 原子 约束 是 同时 的 ， 即 如 
Y= DO= OY X= ASX (2) ”果菜 原子 拟 合 系数 被 约束 为 零 ， 则 所 有 原子 的 该 拟 合 系数 会 被 
式 (2) 中 的 x 为 传 感 算 阵 。 同时 约束 为 零 ， 有 效 的 实现 了 字典 中 的 原子 选择 。 因 此 ， 为 了 
冒号 重 构 内 容 的 关键 是 由 观测 向 量 Y 恢复 @ 。 由 于 多 和 矩 。 使 选择 出 的 原子 更 好 的 实现 信号 的 重 构 , 则 (8) 式 整理 成 目标 : 
阵 已 知 ， 所 以 恢复 9 等 效 为 恢复 原始 信号 X 。 当 M < N 时 ， . 
es ee min|y- AX , +。 (10) 
式 (2) 为 欠 定 方程 ， 所 以 信号 的 重建 为 求解 该 方程 了 = DGX 
的 最 小 六 范 数 解 : 2.2 ”模型 求解 
min|el stimO=7 (3) 式 (10) 是 一 个 基于 2 范 数 的 优化 问题 。 首先 将 其 转换 为 
式 (3) 是 NP 问题 ， 为 了 降低 复杂 度 ， 可 以 用 1 范 数 代 蔡 。 二 阶 锥 规划 (SOCP) 问题 ， 采 用 内 点 法 进行 求解 ， 但 该 方法 计 
minlel szoe@=7 (4) ” 算 复杂 ; 2009 年 ， 文 献 [21] 采 用 梯度 下 降 方法 进行 求解 ， 但 是 
而 OMP 方法 是 压缩 感知 中 经 典 的 贪 楚 迭 代 算 法 之 一 ， 贪 “只 能 解决 特定 问题 。 利 用 文献 0 的 方法 进行 求解 。 
禁 和 迭代 算 法 不 能 直接 解决 最 优化 问题 ， 只 能 找到 近似 最 优 解 。 首先 ， 转 换 式 (10) 为 其 等 价 形式 : 
以 后 的 改进 算法 都 是 在 此 基础 上 优化 而 成 。 了 
min 一 (11) 
2 基于 4 范 数 选择 的 图 像 分 块 稀疏 重 构 ， 
邻 E= 一 |y 一 AX| ， 则 (11) 式 转换 为 : 
2.1 基于 范 数 约束 的 原子 选择 方法 “= ” 
剩余 量 尺 =y， 完 备 字典 A=[a,ay,.…,aw|，3 
假定 初始 化 剩余 量 R=y 字典 A=[@,0,…aw]， 并 miplel, +lxl, siAx+rE=y (12) 
且 过 完备 字典 中 A 的 每 个 原子 为 4 范 数 单位 向 量 即 |a|=1, 根 | 
据 匹 配 追 踪 算法 的 思想 ， 在 过 完备 字典 4 中 选择 出 最 匹配 的 原 整理 后 得 
子 列 w ， 反 复 选 代 以 使 剩余 量 Ro 达到 最 小 。 | ， 由 
min st X7 7 |=y (13) 
观测 向 量 》 可 以 表示 为 人 El 下 
Ns 单位 矩阵 
y= (Ra)a+R, (5) 其 中 / 为 和 位 身 并 _ 
Jj=1 A 
令 B=|X7 rT|,U= 以 式 (13 ) 可 以 优化 成 下 列 形 
RR, 再 按 上 式 进行 分 解 ,直到 余 量 R, 接 近 0 led | ee 
通过 求解 最 小 二 乘 解 可 即 得 原 信号 的 近似 解 式 : 
t=argmin|y-AX|, (9 min|wl ,stBU => (14) 
为 了 减少 重 构 误差 , 从 完备 字典 中 选 出 的 原子 再 用 4 范 数 写 出 式 (14) 的 拉 格 朗 日 算 子 : 
09 进 行 选择 。 
L(VU)=|U|, -7r(A’ (BU -— 15 
a (0)=|wl, -Tr(A (BU -»)) (15) 
i 对 式 (15) 进行 求 导 ， 并 令 其 导数 等 于 零 得 : 
|wl， 二 之 号 = 之 | (7) aL{(U 
lO) pl (10) 
其 中 :矩阵 M = (mh) i 代表 行 号 ，j 代表 列 号 ， 分别 用 mm，m oU 
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信号 X 的 稀疏 度 K 
乘 以 BD- ， 又 因为 BU =y ， 故 ， 输出 : 重 构 信号 元 


1 
其 中 呈 是 一 个 对 角 和 矩阵 冲 ， 对 角 元 素 取 值 为 ““ 下 中 对 式 


(16) 等 号 两 边 


We 


a a) 初始 化 。 将 大 小 为 1, x1. 的 原始 图 像 分 割 为 互 不 重 半 的 
A=2(BD'B’) y (17) ， 
BxN 的 子 块 x (p,q9),p=1,...,B,q 1...,NV, 用 “过 表示 其 中 
结合 式 (17) 推出 ; i,s=1,2.…,S,S =(7, x1,)/(BxN)in=y,A=®@,J=@, 支撑 集 
U=D'B8"(BD'B’) y (18) Tk 
和 迭代 次 数 +=1。 
根据 式 (14) 是 凸 问题 , 所 以 求解 该 式 得 到 的 U 是 全 局 最 优 b) 计算 相关 系数 轧 
解 。 根 据 式 (18) 可 以 看 出 U 与 DD 的 取 值 是 相互 影响 的 。 文 献 c) 第 一 次 第 选 原子 。 uw 进行 降序 排列 , 将 前 大 个 最 大 值 的 
[20] 提 出 了 一 种 迭代 求解 算法 ， 具 体 迭 代 求 解 过 程 如 下 : 索引 存 入 J 中。 
输入 : BeR”™”，yeR”” d) 第 二 次 筛选 原子 (2 范 数 选择 原则 筛选 ) 。 对 7 中 索引 
步骤 : 设 1=0， 初 始 化 D 为 单位 矩阵，D eR 值 对 应 的 原子 的 相关 系数 进行 2; 范 数 选择 ， 选 择 出 对 应 的 
重复 如 下 操作 ， 直 至 收敛 : 的 索引 值 存 入 4 。 
计算 矩阵 Ui,i ， e) 更 新 支撑 集 4: 其 中 4=4DwA4。 


全 信号 逼近 : x=argmin|y 一 Dx 。 


Un=D,'B" (BD,'B’) y 
g) 残 差 的 更 新 : 1+1。 

计算 对 角 和 矩阵 口 ,, ， 其 中 第 i 个 对 角 元 素 为 : b) 判断 远 代 停止 条 件 :每 进行 一 次 远 代 操作 ，++1。 如 果 

i 1<k ， 则 返回 步 又 b) ， 否 则 执行 步 又 i) 。 

~ 21h i) 获取 下 一 个 子 图 像 块 忆 , ， 将 数据 重新 初始 化 ， 

0 太 =》A= 中 ,名 = 中 4J= 中 , 迭代 次 数 !=1， 重 复 步骤 b hh) 。 
we 让 整合 重 构 后 的 子 图 像 块 %,i=1…,5, 恢复 图 像 块 为 重 
构图 像 x 。 


A 
求 出 U 后 ,根据 U -4 便 可 求 出 A， 其 非 零 行 对 应 的 便 


3 ”实验 结果 分 析 


是 贡献 率 较 大 的 原子 列 ， 同 时 也 进一步 有 效 剔 除 无 用 的 干扰 原 
可 本 节 对 所 提 基 于 2 范 数 原子 选择 的 图 像 分 块 稀疏 重 构 算 
2.3 ”基于 4 范 数 选择 的 图 像 分 块 正 交 匹配 追踪 法 进行 实验 仿真 。 所 有 实验 均 在 MATLAB2013a 实现 ， 运 行 平 
将 4 范 数 选择 应 用 到 压缩 感知 中 的 原子 选择 , 结合 分 块 的 台 为 Windows 7 Inter Core i 2.8 Hz CPU 4GRAM。 实 验 采 用 
思想 实现 对 图 像 的 重 构 ， 本 文 称 之 为 641 一 8_OMP 算法。 Lena, Barbara, Camera man, boat 图 ， 尺寸 大 小 均 为 256x256。 
对 图 像 进行 分 块 ， 解 决 了 现 有 算 子 测量 矩阵 存储 量 大 的 问 测试 图 像 选择 Lena 图 ， 通 过 与 匹配 追踪 系列 算法 中 的 四 
题 ， 又 兼备 ROMP 算 子 收敛 速度 快 的 优点 。 种 算法 (BP，OMP，CoSaMP，StOMP) 进行 对 比 。 
本 文 所 提 方 法 以 图 像 块 为 单位 进行 处 理 ， 较 处 理 整 幅 图 像 所 提 方 法 中 ， 采 用 离散 小 波 变换 的 稀疏 化 方法 。 实 验 中 ， 
相 比 ， 容 易 实 现 图 像 重 构 。 采用 高 斯 随机 和 矩阵 为 测量 矩阵。 本 文采 用 分 块 的 思想 结合 范 
DB_OMP 算法 保留 了 ROMP 中 的 正则 化 原子 选择 标 数 思 想 对 图 像 进行 重 构 。 
准 , 并 且 采 用 2 范 数 选择 原则 对 支撑 集 进 行 二 次 筛选 ,具体 的 ， 当 采 样 率 M/N=0.5 时 ， 重 构 算 法 (BP, OMP，CoSaMP， 


次 迭代 过 程 ， 进 行 两 次 原子 筛选 。 第 一 次 原子 筛选 ， 引 入 能 。 StOMP，L 1 一 8_OMP ) 对 Lena 图 像 在 同一 采样 率 下 的 重 构 效 
思想 ， 选 择 多 个 原子 作为 候选 集 。 第 二 次 原子 筛选 ， 从 候选 。 果 对 比如 图 1、2 所 示 。 


证 


每 

集 里 按 2 范 数 选择 原则 考虑 原子 间 的 相关 性 , 剔除 了 干扰 原子 。 由 图 1 可 以 直观 地 看 出 4 一 8_OMP 算法 的 重 构 效果 明显 
快速 、 高 效 选择 出 了 代表 性 原子 并 入 最 优 的 支撑 集 。 优 于 其 他 四 种 重 构 算法 。 从 直观 角度 方面 看 ， 与 图 (b) ~ 〈e) 

以 图 像 块 为 处 理 单位 , 采用 2. 范 数 选择 可 以 有 效 地 剔除 图 。” 相 比 ， 图 (f) 中 帽子 和 头发 的 细节 很 好 地 重 构 出 来 。 

像 块 内 无 用 原子 ， 增 强 快 特征 的 判别 性 。 且 兼顾 了 图 像 块 间 的 2 中 显示 的 是 信号 精确 重建 概率 与 信号 稀 玻 度 玉 之 间 的 
相关 关系 ， 即 某 个 原子 在 某 一 图 像 块 判 定 是 无 用 的 ， 则 该 原子 ”关系 。 在 这 个 实验 中 ， 选 取 重 建 信号 为 可 稀疏 信号 ， 其 长 度 为 
在 所 有 的 图 像 块 上 都 是 无 用 的 ， 回 避 了 图 像 块 间 特 征 选 择 的 不 256; 测量 信号 的 长 度 为 128。 选 择 稀 疏 度 在 10~70 进行 仿真 。 
稳定 性 ， 提 高 了 重 构 的 准确 率 和 效率 。 从 图 2 中 可 以 看 出 ，2. 一 2 -OA 算法 在 稀 疏 度 大 于 50 时 , 仍 


4,1-B_OMP 算法 的 步 又 如 下 : 可 以 工作 ， 其 他 的 算法 作用 甚 微 。 
输入 : 测量 矩阵 BeR””，; 
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选取 Lena(256x256) 


图 片 进行 实验 。 表 1 显示 的 是 当 采 样 


率 相 同时 ， 本 文 所 提 方 法 与 其 他 几 种 算法 在 两 个 


体 指 标 ( 运 


行 时 间 、 峰 值 信 噪 比 (PSNR)) 上 的 对 比 情况 。 运 行 时 间 方 面 ， 
ROMP 所 需 时 间 最 短 ，CoSaMP 所 需 时 间 最 长 。 因 为 CoSaMP 


加 入 下 


溯 的 思想 ， 耗 时 较 长 ，41 一 8_OMP 算法 在 运行 时 间 上 


处 于 中 等 水 平 ， 在 峰值 信 品 比方 面 ，L21 -8_OMP 算法 最 高 ， 
中 只 需要 保存 


大 幅度 领先 其 他 算法 。 这 是 因 


为 在 分 块 月 


E 缩 感 入 


一 个 n*(B*N) 的 矩阵， 而 不 是 一 个 大 和 矩 阵 。 显 然 ， 当 B*N 比较 


小 的 时 候 ， 内 存 


的 干 


用 少 且 计算 速度 快 ; 当 B*N 比较 大 的 时 候 ， 
图 像 的 重 构 效果 比较 好 。 并且 采用 范 数 原子 选择 剔除 支撑 集 
扰 原子 ， 保 证 图 像 块 之 间 相 同 无 用 原子 最 大 程度 的 消除 ， 


提高 图 像 重 构 质量 。 综述 所 述 ，h1 一 8_OMP 算法 融合 2 范 数 


分 块 思想 ， 重 棒 性 强 ， 性 
表 1 采样 率 为 


FE 能 佳 ， 优 势 明显 。 
对， 各 算法 的 的 性 能 


50% 


对 比 
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对 压缩 图 像 进行 分 块 习 
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基于 12,1 范 数 


E 构 时 ， 先 将 实验 图 像 分 成 B*N 个 小 块 ， 
分 块 示例 如 图 3、4 所 示 , 并 且 保 证 对 每 个 小 块 采 取 相 同 实 验 处 
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(c)4*2 块 
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图 3 Lena 在 不 同 分 块 下 的 重 构 效 果 


(a)2*1 块 


(e)8*1 块 


令 对 应 的 讽 
后 的 向 量 表示 为 


i 


从 上 式 可 


n*(B*N) 的 矩 阵 ， 而 不 是 一 个 大 和 矩 阵 。 


(b)2*2 块 


(8*2 块 


时 候 ， 内 存 


5 


j 少 且 计 算 速 度 快 ; 当 


y i= Pax 
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图 4 Barbara 在 不 同 分 块 下 的 重 构 效果 


(d)4*4 块 


(hb8*8 块 


量 值 为 Y， 则 原始 图 像 的 第 i 个 小 块 经 过 变换 


以 看 到 ， 在 分 块 压缩 感知 中 只 需要 保存 一 个 
显然 ， 当 B*N 比较 小 的 
B*N 比较 大 的 时 候 ， 图 像 


的 重 构 效果 比较 好 。 所 以 需要 找到 一 个 比较 折 中 的 B*N 的 值 。 


然后 , 对 每 个 小 块 采 | 


基本 


实现 图 像 重 构 。1 
的 块 效应 明显 。 


4 ”结束 语 


本 文 提出 了 一 种 结合 2 范 数 选择 


运行 时 间 /s 


5.0034 2.2554 5.2234 2.6754 8.2100 2.4432 


峰值 信 噪 比 /db 26.4102 18.9943 27.5643 25.5404 26.7988 25.3743 


5.2446 


67.3043 


配 追 踪 的 重 构 方法 。 该 算法 以 块 结构 获 ] 


Fb 范 数 选 择 的 正 交 匹 配 追 踪 ， 从 而 


图 3、4 看 出 (a)(b)(c)(e) 视 觉 效果 良好 ， 剩 余 


表 2 可 知 ， 块 状 数量 选取 不 同 ， 习 
可 知 ， 在 分 块 数量 为 4*2 时 ， 人 恢复 效果 佳 。 


| 四 
a 
Ey 
Ne 
洗 
mm 
湘 
习 


同 ， 由 统计 


原子 的 基于 分 块 正 交 匹 
取 图 像 ， 利 用 


4 范 数 最 


小 化 的 选择 方法 选择 出 各 个 分 块 中 的 代表 原子 ， 有 效 地 剔除 了 


为 了 方便 图 像 进行 传输 ， 本 文采 取 对 


图 像 进行 分 块 操作 。 


干扰 原子 ， 降 低 计算 量 并 且 
测试 条 件 一 致 的 前 提 下 ， 从 主观 视觉 和 客观 数据 两 方面 


构 效 果 好 。 实 验 结果 表明 ， 保 证 


i 都 说 明 


录用 稿 


了 该 算法 性 能 的 有 效 性 。 
其 


| | 《 V 加 :由 a 
朱 华 ， 等; 范 数 原子 选择 ; 


然而 ， 当 受 污染 的 数据 分 布 在 块 状 结 


构 的 边缘 处 时 ， 其 重 构 效 果 有 所 下 降 ， 故 在 这 一 方面 还 有 待 提 


出 改进 措施 。 
表 2 采样 率 0.5， 不 同 分 块 下 的 的 性 能 对 比 
Lena Barbara camera man boat 

Psnr Time/s Psnr Time/s Psnr Time/s Psnr Time/s 
2*] 66.7480 6.0744 26.8162 79.3814 36.4170 19.2552 47.3569 76.3372 
2*2 63.6056 9.4125 27.2386 71.4876 37.1783 18.1028 44.2096 71.1920 
4*2 67.3043 5.2446 32.1075 71.6157 40.4319 20.0217 46.7853 71.1635 
4*4 65.0593 4.1227 28.2386 132.7240 39.1449 28.7785 44.6225 115.2722 
8*] 65.3982 9.6814 30.4458 202.5869 36.3285 32.9373 39.3478 138.8600 
+ 人 2 13.0529 9.5404 8.2571 244.6040 11,2366 34.1345 16.7334 133.9920 
8*4 12.7331 6.1774 7.3321 219.0565 10.5632 20.3923 11.5433 173.4998 
8*8 11.2948 5.6814 4.5861 179.8573 9.2342 37.4083 9.2234 193.3932 
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